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We calculate the parallel (VV) and perpendicular (VH) polarized Raman spectra of amorphous 
silica. Model Si02 glasses, uncompressed and compressed, were generated by a combination of 
C , classical and ah initio molecular-dynamics simulations and their dynamical matrices were computed 

C . within the framework of the density functional theory. The Raman scattering intensities were 

' determined using the bond-polarizability model and a good agreement with experimental spectra 

was found. We confirm that the modes associated to the fourfold and threefold rings produce most of 
the Raman intensity of the Di and D2 peaks, respectively, in the VV Raman spectra. Modifications 
of the Raman spectra upon compression are found to be in agreement with experimental data. We 
show that the modes associated to the fourfold rings still exist upon compression but do not produce 
a strong Raman intensity, whereas the ones associated to the threefold rings do. This result strongly 
suggests that the area under the Di and D2 peaks is not directly proportional to the concentration 
of small rings in amorphous Si02. 
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PACS numbers: 61.43.Fs,63.50.+x,78.30.-j,71.15.Pd 



INTRODUCTION 



> 

Raman spectroscopy has provided a powerful tool to study the dynamical properties of amorphous materials. 
00 ■ The interpretation of the vibrational spectra of vitreous silica (a-Si02) has been of interest for many years and 
' considerable progress has been made. Despite this fact, the origin of the Di and D2 lines at 495 cm^^ and 606 cm~^, 
: respectively, in the Raman spectrum is still not fully understood. In order to explain the frequencies of these bands, 
Galeener et alm'^'^'^ has argued that D2 can be assigned to the breathing modes of a planar threefold ring, and Di to 
that of a regular but slightly puckered fourfold ring within the continuous-random-network (CRN)iSii. On the other 
hand and within the submicrocrystallite approach, PhillipsSiSiifliii has proposed that within the a-Si02 structure 
C I there exist two types of non-bridging oxygen (NBO) sites at the surface of the cluster, which can be associated to the 
1-^ ■ structure responsible for the Di and D2 defect lines. To date, much work has been done in the light of both models. 
P5 ' Using a pair potential and the bond-polarizability (BP) mode&i^, Murray and Chingi^ have calculated the Raman 
O • spectra for normal and compressed Si02 glass models based on the CRN model with periodic boundary conditions. 

' Their results generally reproduce the shapes of the bands present in experimental Raman spectra, particularly the 
^ : Di and D2 lines in contrast with the paracrystalline theory of Si02. Moreover, the ring statistic analysis showed 
that there are no three or fourfold rings in their models. This means that the Di and D2 lines are not strictly due 
. _ to the threefold rings and fourfold rings as suggested by Galeener. To conciliate this later point of view and their 
■ work, Murray and Ching indicated that the Si-O-Si angles close to those found in the three and fourfold rings are 
" ■ ' responsible for these modes. More recently, using molecular-dynamics (MD) simulations and the BP model, Zotov 
et alm^ have calculated the Raman spectra of a-Si02 models. The obtained VV Raman spectra were, generally, 
in good agreement with experiment but, despite the presence of three and fourfold rings, their models failed to 
reproduce the Di and D2 peaks. Using a first-principles density functional approach, Umari et aii^ have calculated 
the Raman intensities of a-quartz by evaluating the variation of the polarizability tensors for finite displacements of 
the atoms. This work supports the use of the bond polarizability model for the calculation of the Raman intensities 
in a-Si02. By applying this method to an amorphous Si02 sample containing a large concentration of threefold 
rings, they were able to relate the intensity of the Di and D2 peaks to the concentration of four and threefold 
rings, respectivelyiSi. In a recent study, Benoit and Kob compared the vibrational properties of model Si02 glasses 
generated by molecular-dynamics simulations using the effective force field by van Beest et al (BKS)iI with the ones 
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obtained if the BKS structure is relaxed using ab initio calculationslS.. It was found that this relaxation significantly 
improves the agreement of the density of states with the experimental result. Therefore it has become of interest to 
use the data of this approach to determine the Raman spectra of these systems. 

In this work, we present the polarized Raman spectra of compressed and uncompressed a-Si02 samples prepared 
by combined classical and ab initio molecular dynamics simulations. Raman spectra were calculated using the bond- 
polarizability model and diagonalisation of the dynamical matrix, as obtained from the evaluation of the second 
derivatives of the total energy with respect to atomic displacements by taking finite differences of the atomic forces. 
In the following section, we describe the models used to compute the Raman spectra. Then in the third section we 
present the calculated spectra and a detailed analysis of the Raman signatures of the small rings. Finally the spectra 
are also calculated for compressed silica and a discussion in terms of the Raman signature of the rings is also given 
before the conclusion. 



II. COMPUTATIONAL METHOD 



The time averaged power flux of Raman scattered light in a given direction, with a frequency between tof and 
ujf + dujf in a solid angle dfl, is related to the differential scattering cross-section: 



+ l)h ^ VaVpHajp\{uj)WjWx, (1) 
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where lu = Ui ~ luj and 
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a„,(j) = 5]-^<n^|j>, (3) 

LOi is the frequency of the incident light, v and w are the polarization unit vectors for scattered and incident light 
respectively, B{lo) is the Bose factor, Af„ the mass of the nth atom, ujj and < n6 \ j > axe respectively the frequency 
and the {n5) component of the jth mode. The Greek characters denote the Cartesian components (x,y,z). The 
coefficients tt^^ ^ connect the polarization fiuctuations to the atomic motions^^'^^ and they are obtained by expanding 
the polarizability tensor rf" in terms of atomic displacements u^, with 
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In addition to these coefficients, the calculation of Raman intensities requires mode frequencies and eigenvectors: Uj 
and < I j >, obtained by solving the equation: 

^|j>=c^f|j>, (5) 

where D is the dynamical matrix of the system, with the element Da/jinm) —< na \ D \ to/3 >. Each element of the 
dynamical matrix, D, is given by: 

Dapin^m)^ I (t>ap{n,m), (6) 

with (pap^UjUi) being the force constants between atoms n and to. In our calculations, the interactions are given by 
a first-principles approach. 

If we assume that scattering can be described within the framework of the bond polarizability (BP) theory, the 
polarization is only modulated by nearest-neighbor bonds and the components of the induced polarizability tensor tt 
are given hyi^ 

nap{r) = ^{ai + 2ap)Sap + {ai - ap){faff3 - ^S^p), (7) 
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where f is the unit vector along vector r which connects the n and m atoms linked by the bond. The parameters ai 
and Up correspond to the longitudinal and perpendicular bond polarizability, respectively. Within this approach and 
in order to evaluate the derivatives, one assumes that the bond polarizability parameters are functions of the bond 
lengths r only. The derivatives tt^^ ^ are then given bjsi^ 

^2/3,7 = X! "^(^"P + "i)^Q/3^7 + ("i - a'p){faf0 - \5al3)f.y + — ^(^a^r^ + Sp^f^ - 2faf0f^) (8) 

m 

where a' — {da/ dr)\^^^^, vq is the equilibrium distance. 

In our a-Si02 model, a single type of bond occurs (Si-0) and the bond polarizability model is completely defined by 
three parameters: a = {2ap + a[) , j3 = a'l — a'p and 7 = (a/ — ap)/r. 



III. RAMAN SCATTERING IN AMORPHOUS SIO2 

The amorphous Si02 samples used in this study contain 26 Si02 units in a cubic box with periodic boundary 
conditions at the experimental density of 2.2 g.cm"'^. They were generated by molecular-dynamics simulations using 
the van Beest et al. force fieldlii: Three well-equilibrated liquids at 3000 K were quenched to 300 K using three 
different quench rates (7 • 10^° K/s (A), 3 • 10" K/s (B) and 5 • lO^^ j^/s (C)). The structures of the model samples 
were then refined using first-principles calculations of the Car-Parrinello type at 300 K using the CPMD code^^ (for 
the simulation details, see Rcf.^*). After a relaxation of the atomic positions to K in the density functional theory 
framework, the dynamical matrices of these samples were computed by evaluating the second derivatives of the ab 
initio total energy with respect to the atomic displacements by taking the finite differences of the atomic forces. The 
vibrational frequencies and corresponding modes were obtained by diagonalisation of the dynamical matrices. 

The structural characteristics of the three samples were first analyzed in terms of the pair correlation functions 
and of the static structure factors. All samples present a perfect tetrahedrally connected network and no noticeable 
sample specificities were found on such averaged quantities. 



A. Density of states 

The vibrational properties of a-Si02 have been extensively investigated by many experimental and some theoretical 
means. In order to show the validity of the vibrational features of our model, we report, in Figure ^ the calculated 
effective neutron scattering cross section G{u!) = C{uj)g{uj), where g{uj) is the vibrational density of states and C{uj) 
was computed using the incoherent approximation, following Ref.^^. The calculated G{uj) was obtained by averaging 
over the G{uj) of the three Si02 samples A, B and C. The calculated G{uj) is compared to the experimental one— 
and a good agreement between the two is observed. Shapes and positions of the principal peaks are well reproduced 
with a low-frequency band at 420 cm^^, an intermediate frequency band at 790 cm~^ and a double peak at 1050 
cm~^ and 1195 cm~^ in the high-frequency regioniSii^iii. Note that the lack of a small peak at around 4 THz in the 
experimental data is related to the insufficient experimental resolution—. 



B. Raman spectra 

The purpose of the present subsection is to present our calculation results for the polarized parallel (VV) and 
perpendicular (VH) Raman spectra of different samples of a-Si02. We assumed that scattering was produced 
by the BP mechanism. In all spectra, intensities represent the so-called reduced intensities, which are obtained 
by multiplying the experimentally measured spectrum with the correction factor (lu/{B{lu) + 1)). In the VV 
configuration, both the incident and scattered polarizations are along the z axis and, for the VH configuration, the 
incident and scattered polarizations are along the z and x axis respectively. The intensities of the different spectra 
are normalized and, to make the comparison with the experimental results more realistic, we averaged over all sample 
orientations in the Raman spectra calculations. 

Firstly the Raman spectra were calculated for normal density (2.20 g/cm^) a-Si02 samples. The empirical param- 
eters used in this study are those derived from the results of Umari et al.^^ for a-quartz: q;=1.0, /3=0.25 and 7=0.07. 
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The calculated Raman spectra (solid lines), which are obtained by an average over the three different a-Si02 samples, 
are plotted together with the experimental result of Galeener— in Figure |21 (dashed lines). 

The calculated VV Raman spectrum (FigureEl^a)) is dominated by a lower band, around 500 cm^^ and reproduces 
the Di — 495 cm~^ and D2 = 606 cm~^ peaks found in experimental spectraS4ii*2SiS&. The calculated Di and 
D2 occur at 489 cm^^ and 614 cm~^, respectively. For these bands, the curve is very close in peak positions and 
qualitatively in peak shapes to the experimental curve. We note also that the bands at 436 cm^^ and 790 cm^^ 
are close to the experimental bands observed around 430 cm~^ and 800 cm~^iSiii*2£iSSi21. In our calculation within 
the BP approximation, the double peak found experimentally in the high frequency region (1100-1200 cm~^) are 
better reproduced than in previous worksi^ii^. Concerning the VH spectrum. Figure |2b), three major bands can 
be distinguished. The double peak in the high frequency region (1100-1200 cm~^) are well placed and the curve is, 
generally, close to the experiment. The lower band (< 500 cm~^) and the intermediate one (around 793 cm~^) are 
positioned correctly but their intensities are too weak in comparison with the experimental lines. 

In the VV Raman spectrum, the peaks around 430 cm^^, 800 cm^^ and 1100-1200 cm^^ have been understood 
in terms of the vibrations of a CRN modelli^. However the nature of the "defects" responsible for the Di and 
D2 bands in the VV spectrum of a-Si02 is still the subject of a theoretical controversj^'- . As we have recalled 
above, according to Galeeneji^^iMi^j the Di peak is due to the vibrational modes associated with oxygen atoms 
belonging to planar fourfold rings and the D2 peak to the ones associated with oxygens belonging to threefold 
rings within the CRN modcl^-^. In contrast to this, PhillipgSiSiiSiii argued that these bands may be associated 
with skeletal surface Si02 vibrational modes in the framework of paracrystalline theory. In the present work, the 
structural characteristics of the studied systems showed, on one hand, that they contain threefold and fourfold 
rings and on an another hand that, as the models of Murray and Chingi^, none of them contain NBO sites or 
free surfaces. In Figure |21 we present the different ring size distributions of the three amorphous samples A, B 
and C. As it can be seen, the shapes of the ring size distributions are quite different for the three samples. Since 
no particular trend can be derived with respect to the quench rates used to generate these model structures, 
the differences found in the distributions arc attributed to the statistics. However it is interesting to note that 
there is a significant amount of three and fourfold rings in all three samples. In previous molecular-dynamics 
simulations^'*, it has been shown that the concentration of small rings increases with an increasing quench rate. As 
a consequence, the concentration of small rings in our samples is certainly much larger than in real Si02 samples. 
We find that the number of threefold rings varies from 1 in samples A and C, to 2 in sample B, and the number of 
fourfold rings are equal to 4 in sample A, to 3 in sample B and to 6 in sample C. Nevertheless, as shown in Figure 
13 the VV Raman spectra of our models present both Di and D2 peaks which means that the later may be cor- 
related with the presence of the three and fourfold rings as argued by Galeener and excludes the Phillips point of view. 



C. Raman signature of the rings 

In order to clarify whether vibrational modes associated with the three and fourfold rings can be responsible for the 
Di and D2 lines, we analyzed the Raman intensity due to these specific geometries. However the connection between 
a given geometry in the amorphous structure and the Raman intensity produced by it, is not straightforward at all. 

In this perspective, we carried out an analysis on the three different samples individually. We first analyzed the 
modes that gave the strongest Raman activity in the 470-515 cm"-'^ range for the Di line and in the 550-730 cm~^ 
range for the D2 line. The selection of these modes was done first by ordering the modes as a function of their Raman 
activity in the chosen frequency region. We then added together the activities of the most Raman-active modes until 
the sum reaches 80 % of the total sum: Only the modes included in those 80% were subsequently studied. 

Then, for the selected modes, we computed the localization entropy defined as: 

no 

where \j > are the eigenmode components of mode j and \nd > are the S eigenmode components of atom n. The 
above defined entropy is connected with the lack of information about the position of the jth phonon and gives the 
same kind of information than the participation ratio. If the mode is perfectly localized, then the lack of information 
is minimum and only one of the | < nS\j > ^ is equal to one, which leads to Sj = 0. On the other hand, if the jth 
mode is completely delocalized, then | < nS\j > p = 1/3A^, V nS and Sj is maximum, i.e. equal to 1. Here Sj has 
been normalized to 1 for convenience, thus it is not a genuine entropy. Let us consider the case where the jth mode 
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is localized on a given ring composed of Nji atoms. The value of Sj would then be around: 

nSRing 

Thus if a mode is locaHzed on all atoms in a threefold or a fourfold ring, the corresponding Sf- would be close to 0.529 
or 0.582, respectively. Table Ogives the localization entropies of the most Raman-active modes for every sample in 
the chosen frequency ranges. Note that Sj gives an estimation of the localization of a mode but does not tell whether 
a mode is localized on a specific ring. 

From tabled it is clear that the most active Raman modes in the 470-515 cm~^ frequency range do not correspond 
to very localized modes since their localization entropies lie around 0.8 and are never smaller than 0.67. The modes 
responsible of the Di line can therefore not be considered as very localized modes. For D2, the situation is similar 
since the localization entropies lie also between 0.72 and 0.85 for the most intensive modes in the 550-730 cm^^ 
frequency range. Therefore this critcrium does not appear to be useful in order to understand the nature of the 
modes responsible for the Di and D2 lines. 

We then tried to find other criteria that could help to determine whether these modes can be associated to the 
three and fourfold rings or not. 

The first chosen criterium is based on vibrational considerations: For every mode and for each ring, we computed 
the following quantity 

ARi„g(j)- E I (11) 

nGRing 

If the atomic displacements due to the considered mode are of similar amplitude for every atom, then the value of 
^Ring(j) would be equal to 3 x 8 / 3A^ for the fourfold rings and to 3 x 6 / 3A^ for the threefold rings since there 
are 8 or 6 atoms included, respectively. We therefore computed the quantity i^Ring(j) = ^Ring (j) X iN/2,Nii where 
iV/j=8 for the fourfold rings and Nji=6 for the threefold rings, in order to obtain an estimation of the participation 
of atoms belonging to a given ring, for a given mode. If the jth mode is completely localized on a ring, one obtains 
the maximum value for Pring(j): -fpungO) — ^•'''^ fourfold rings and P^i^gij) — threefold rings. Let 

us recall however that for a completely delocalized mode, PmngU) is equal to one and that Pmngij) = means that 
the considered ring is not concerned by the jth mode. 

The second chosen criterium is based on the contribution of a given mode and a given ring to the total Raman 
spectrum. This contribution is evaluated by computing what we name the "Raman ratio" which is defined in order 
to estimate the contribution of each ring to the Raman intensity of a given mode. Let us consider a given ring. The 
Raman intensity can be divided in three contributions by considering separately the atoms belonging to the ring and 
the other ones, in the following way: 

.„,(j) = ^<nS\j>+ ^<nS\j> (12) 




(13) 

For a given polarization (we omit the Greek indices), the contribution to the Raman scattering can be written as: 
Htot^uj) = ^a*^(j>/3A(j)^(^(t^-t^j)-i5(w-hu;j)) (14) 

= E K7(J) + <M* Ka(j) + ^pxij)] ^(^(^ - - + .;,)) (15) 

= E KijHxU) + c,(j>;^a(j) + {«",(j>;;a(j-) + «";(j>^a(j)}] (i6) 

x^('5(^-^j)-'^(w + ^j)) (17) 

= Hnot Ring(t^) + -ffRing(t^) + -^overlap (t^) (18) 

To obtain the contribution of the atoms belonging to a given ring to a selected mode jo, we computed the integrated 
intensity of the spectrum, averaged over all possible sample orientations, in the chosen frequency range: 

^ing — ^otal — V^°^. j^jjjg (19) 
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with 



and 



Vtotai = J < ■f^totai(w) > duj (20) 



not Ring 



/ < ^not Ring(^) > (21) 



where the < ... > denote the average over all possible orientations and Ringl'^) been obtained by setting to 
zero the < nd\jo > atomic components of the considered ring for the selected mode jo- The integrals run over the 
frequency ranges corresponding to the Di line or to the D2 line, depending on which type of rings we want to consider 
(fourfold or threefold). 

We then needed a quantity that could be used in order to compare the contributions of different rings to the same 
mode. For this, we defined the relative Raman participation ratio as follows. The integrated Raman intensity of the 
spectrum, averaged over all possible orientations, has been computed without the contribution of mode jo^ 



where 



and 




The relative Raman participation ratio is then obtained by computing the quantity: 

Raman ratio(jo) = (25) 

The Raman ratio gives the relative participation of the atoms belonging to a given ring to the Raman intensity 
produced by the selected mode jo- 

If the atoms of the chosen ring do not participate to the mode jo then V^°^^ will be equal to zero and the Raman 

ratio as well. On the other hand, if only the atoms of the ring participate to the Raman scattering, then V^°^^ = T4otai 
and the Raman ratio is maximum, i.e. equal to Vtotai /(^4otai — V^ot jo)- However, from the definition of V^^jg (Eq. 
I19|l . one can see that the Raman ratio can be negative: This is due to the fact that the -ffoveriap term can be negative 
with a magnitude larger than that of -f^Ring for the chosen ring. A negative Vp°„g would mean that there is a strong 
negative "coupling" between the atoms of the ring and the other atoms. A similar approach has already been used 
by Lazzeri et al. for the study of the ring contributions to the Raman spectra of crystalline silica phases and zeolites 
and has produced negative overlaps in these phasesSS. On the other hand, a large and positive Raman ratio would 
mean that the atoms of the ring participate to the mode but also that the coupling with the other atoms might be 
large and positive. 



1. Di line 

The results for the most Raman-active modes in the 470-515 cm~^ frequency range are presented in table ITll for the 
A, B and C samples with the vibrational amplitudes (PRing) on the left columns and the Raman ratio on the right 
columns. 

If -PRing is larger than one for a given ring, the atoms belonging to this ring have larger displacements than the 
other atoms for the considered mode. From table |nl one can clearly see that the selected mode always gives rise to 
larger displacements on atoms belonging to one (or more than one) of the fourfold rings. As a consequence, one is 
tempted to "associate" Raman-active modes to specific rings: For instance, the 474 cm^^ mode in sample B generates 
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large displacements on atoms belonging to ring 1 in that sample and the 507 cm"-'^ mode in sample A generates large 
displacements on atoms belonging to rings 1 and 2 in that sample. 

The fact that some modes can be associated to several rings lead us to check whether these rings had common 
atoms. The result of this analysis is presented in tabic UTTl Most of the time, when modes give large displacements 
on several rings, the later have indeed common atoms, but this is not always true. For instance, the 500 cm~^ mode 
in sample A generates large displacements on rings 2 and 3, which do not have common atoms. This result indicates 
that there exist ring-like modes which concern more than one ring. 

The Raman ratio defined above has been computed for the most Raman-active modes in samples A, B and C and 
for every fourfold ring in these systems. Results are presented on the right sides of tabled] for samples A, B and C 
separately. We first notice that the Raman ratio varies approximately from 1 to -1 and thus can be strongly negative. 
If we consider the largest positive Raman ratio, we note that they almost always correspond to the largest values of 
^'Ring- By cross checking the largest terms in both the vibrational amplitudes (Pmng) and the Raman ratio, we can 
easily deduce that if the Raman ratio is large and positive, the atoms of the ring undergo larger displacements than 
the other atoms and therefore that the contribution of these atoms to the mode is important. On the other hand, if 
we consider the negative Raman ratio, we note that they almost always correspond to small - or close to 1 - Pmng 
values. This result is an indication that if the Raman ratio is negative, then the atoms of the ring do not participate 
much to the mode. 

In order to illustrate the results presented in tables m we show examples of Raman spectra from which the 
contributions of the atoms in individual rings have been extracted. Let us consider the mode at 488 cm^^ of sample 
C. In figure^ the total spectrum {Ivv) and the spectrum where this mode has been removed (Imodc) are depicted 
in bold lines. For this mode, the spectra for which the eigenmode components of atoms in rings 1, 3 and 6 have 
been set to zero, are shown in dashed lines (iRing)- As it can be easily seen, for rings 1 and 6, the Imodc and iRing 
almost exactly coincide (Imodc - iRing is almost zero) and the subtractions lyy-Imode and Ivy-lRing are identical. This 
situation corresponds to a Raman ratio close to one in table^land indicates that the atoms of rings 1 and 6 contribute 
significantly to the 488 cm~^ mode. On the opposite, for ring 3, the iRing gives rise to a larger Di peak than that of 
the total VV spectrum, leading to a negative lyy-lRing contribution and very different Imodc and iRing spectra. This 
indicates that the atoms of ring 3 do not participate to the 488 cm~^ mode. 

As a summary, our results suggest that the modes which produce more than 80 % of the Raman activity in the 
frequency region of the Di line can almost always be associated to the atoms belonging to one or more fourfold 
rings (i.e., the Raman ratio is larger than 0.7). Furthermore more than 70 % of the atomic displacements due to 
these modes (PRing) are located on atoms belonging to the fourfold rings. One exception can be found for the 
496 cm~^ mode in sample C, for which the Raman ratio is not larger than 0.5 for all rings. 



2. D2 line 

Since the D2 line has been suggested to be due to modes involving the threefold rings, we analyzed the Raman modes 
in the corresponding frequency region with respect to atoms belonging to such rings. In the three Si02 samples, the 
threefold rings are less numerous than the fourfold ones (1 in sample A, 2 in sample B and 1 in sample C), hence 
reducing the statistics in that case. 

The results for the most Raman-active modes in the 550-730 cm^^ frequency range (corresponding to the D2 line 
in our spectrum) are presented in t able Hvl for samples A, B and C with the vibrational amplitudes (Pring) on the left 
columns and the Raman ratio on the right columns. As for the Di line, we can observe that there is a correlation 
between the largest Ppung values and the largest positive Raman ratio. We can also note that not all the Raman-active 
modes in the D2 line region correspond to modes associated to the threefold rings, whereas the contrary is true: All 
rings participate to, at least, one Raman-active mode in the D2 line region. 

As for the Di line, we illustrate the results presented in table Hvl bv the Raman spectra depicted in figure [S] In 
this figure, we chose to show the Raman spectrum in which the mode of frequency 609 cm^^ in sample B has been 
removed (Imodc) together with the Raman spectra in which the eigenmode components of atoms belonging to the 
threefold rings have been set to zero (iRing for rings 1 and 2). It is clear that if the eigenmode components of atoms 
in ring 2 are set to zero, the Raman spectrum is identical to the one obtained if the selected mode is removed, i.e. 
the resulting difference iRing - Imodc (dotted line) is zero. On the opposite, the contribution of atoms of ring 1 to the 
609 cm~^ mode is non-existent. 

Recently Umari and Pasquarelloi^ have deduced the proportion of oxygen atoms belonging to the threefold and 
fourfold rings in an experimental sample using a direct relation between the areas under the experimental and cal- 
culated Di and D2 peaks, respectively. In view of the results presented in tables Ull and IIVI such an approach seems 
a bit hazardous. Indeed, we have demonstrated that among the most Raman-active modes in the regions of the Di 
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and D2 lines (i.e. modes that produce more than 80 % of the peaks), some can not be associated to modes involving 
atoms of a fourfold or a threefold ring. This result means that there is no proportionality relation between the area 
under the peaks and the proportion of threefold and fourfold rings in the Si02 glass. 

IV. RAMAN SCATTERING IN COMPRESSED AMORPHOUS SIO2 
A. Generation of the compressed samples 

In order to mimic the effect of pressure, the box size of the amorphous sample C has been rescaled to three different 
sizes in successive steps. First, the box length and coordinates were rescaled from the experimental density of 2.20 
g/cm'^ to a density of 2.40 g/cm^ (sample CI) and a Car-ParrincUo molecular-dynamics simulations of 1.9 ps was 
launched at 300 K, using the CPMD codeSi. The entire procedure was repeated for two other densities: 2.67 g/crn^ 
(sample C2) and 4.28 g/cvo? (sample C3). For the sample at 2.67 g/cm^^, the initial coordinates were taken and 
rescaled from the final ones of the simulation performed on sample CI, and for the sample at 4.28 g/cm^, the initial 
coordinates were taken and rescaled from the final ones of the simulation performed on sample C2. In this latter 
case, an important modification of the structure leading to an large increase of the ionic temperature was observed. 
Therefore during the first ps of the ab initio MD simulation, the ionic temperature has been constrained to stay 
at 300 K. This constraint was then released for the rest of the ah initio MD simulation on sample C3. In order to 
achieve good convergence for the calculation of the pressures in samples CI, C2 and C3, higher energy cutoffs would 
be required. We estimated that this information was not worth the computational effort and obtained the pressures 
from the volume to pressure conversion using the experimental equation of state of a-SiOi'^'^. The pressures of our 
samples are thus evaluated to w 3 GPa for sample CI, « 7 GPa for sample C2 and « 34 GPa for sample C3. 

For samples CI and C2, we observe no modifications of the network topology. The observed structural modifications 
only concern the shift of the Si-O-Si angle distributions towards smaller angles (Figure IHI) and no changes in the Si-0 
bond lengths and in the ring statistics are noticed. In contrast to this, the structure of sample C3 is very different 
from the ones of sample C, CI and C2: The tetrahedral network is destroyed and replaced by a mixing of four, five 
and sixfold-coordinated silicon atoms connected by corners and/or edges, and threefold-coordinated oxygen atoms. 
These coordination changes explain the very different Si-O-Si and 0-Si-O angular distributions observed in Figure 

At the end of these ab initio molecular-dynamics simulations, the atomic coordinates were relaxed to K and the 
dynamical matrices were computed as described in section HTH For every sample, we then obtained the vibrational 
frequencies and the corresponding modes by a diagonalisation of the dynamical matrices. The vibrational densities 
of states of sample C and of the three compressed samples CI, C2 and C3 are depicted in Figure As pressure is 
increased from C to C2, the double peak in the high frequency region disappears and the peaks at ~600 cm~^ and 
~800 cm~^ are shifted to higher wavenumbers. One can also notice a decrease of the peak intensity at ^--^400 cm~^. 
On the other hand, the vibrational density of states of sample C3 exhibits very different features which reflect the 
changes in the structure: No clear gap can be seen in the broad peak going from ~100 cm^^ up to 1300 cm~^. 

B. Raman spectra and signature of the rings 

In FigurelHl the VV (left panel) and VH (right panel) Raman spectra of the compressed samples CI, C2 and C3 are 
compared to that of sample C. Firstly, we note that the VV and VH spectra of sample C3 are very different compared 
to that of samples C, CI and C2, which is due to its very different density of states and structure. For the slightly 
compressed samples (CI and C2) the main band and the Di line in the VV spectra merge into a sharp main peak. The 
intensity of this peak increases with increasing density and moves to higher wavenumbers. The increasing sharpness of 
the band around Di with increasing density could be attributed to the shift to lower angle and the increasing sharpness 
of the Si-O-Si angle distribution shown in FigurelHl In contrast with this behavior, the D2 line intensity does not 
increase with increasing density, however its position slightly shifts towards higher wavenumbers. One can also note 
that all the bands lower than 900 cm~^ clearly shift to higher wavenumbers. On the opposite, in the high frequency 
region, we observe that the double peak slightly shifts to lower wavenumbers. Overall these behaviors have been 
observed experimentalljsSS*^ and reproduced by other numerical studies^^. For the VH spectra, we observe that the 
principal bands shift to higher wavenumbers except for the largest peak of the double peak in the high frequency region. 

In order to show the effect of compression on the ring contributions, we have evaluated the localization entropies, 
the vibrational amplitudes PRing and the Raman ratio as defined in section UlI CI for all rings in samples CI and C2. 
The sample C3 was not studied in that respect since the tetrahedral network is damaged in that sample. The results 
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are presented in table for the entropies and in table I VII and table IVIII for the PRing values and the Raman ratio in 
the frequency regions of the Di and D2 lines, respectively. 

The results presented in table El show that the modes contributing to the sharp line around 500 cm~^ are less 
numerous in the compressed samples CI and C2 than in sample C. For samples CI and C2, the frequency range in 
which these modes have been sought for was modified to take into account the frequency shift of the peaks. These 
frequency regions have been set to 490-540 cm~^ for sample CI and to 505-585 cm^^ for sample C2. It should be 
noted that the intense sharp peak at 540 cm~^ in sample C2 is due to only two modes that produce 80 % of the total 
Raman intensity. However, the entropies of these modes are not much different from the ones of the modes responsible 
of the Di line in the uncompressed samples. A similar conclusion on the localization entropy can be given for what 
concerns the most Raman-active modes in the frequency region around the D2 line (605-770 cm~^). Note however 
that there are more Raman-active modes in this region than in the D2 region of the uncompressed sample C. 

The vibrational amplitudes (PRing) and the Raman ratio (see definitions in section Fill C|) were then evaluated for 
the selected modes in the frequency regions of interest. The results for the modes in the 490-540 cm~^ and 505-585 
cm~^ frequency regions are presented in table lVll The mode at 503 cni^^ in sample CI is clearly the same mode than 
the 488 cm"'^ one in sample C since the PRing values and the Raman ratio are maximum on rings 1 and 6 for this 
mode. The mode at 512 cm~^ in sample CI could also be related to the 496 cm~^ mode in sample C if one considers 
the value of PRing which is maximum on ring 5. Upon compression these two modes have been shifted of roughly 15 
cm~^. However, all the other modes that were associated to one or two rings in sample C, have disappeared in sample 
CI and, upon further compression in sample C2, even the modes that were still present in sample CI vanished. For 
sample C2, the remaining modes in the 505-585 cm~^ frequency range do not present large values of PRing or of the 
Raman ratio: This indicates that the Raman-active modes in this frequency region can not be associated to motions 
of atoms belonging to a fourfold ring. Since the same rings of sample C still exist in sample C2, this result means that 
the compression induced a structural modification of the rings and of their environment that produced the extinction 
of the ring-like modes in the Raman spectra. When sample C is compressed into CI and C2, the network connections 
(and thus the rings) do not change with increasing density but only the angles and relative distances between atoms 
do. The changes in the Raman spectra are therefore not due to a different proportion of rings but only to changes 
in the structure of the rings, mainly the decrease of the Si-O-Si angles. A more detailed analysis showed that the 
ring-like modes still exist in this frequency range (i.e. there are modes for which PRing is maximum for a given fourfold 
ring) but their Raman activities are much smaller, if not zero. 

The effect of compression seems to be a little different for the modes in the frequency region corresponding to the 
D2 line. The PRing values and the Raman ratio for the 605-770 cm~^ frequency range are presented in table IVIII 
for the compressed samples CI and C2. In both samples, one can see that there is always one mode being clearly 
associated to the threefold ring: The 646 cm^^ mode in sample CI and the 664 cm~^ one in sample C2. This result 
indicates that the mode at 619 cm~^ in sample C is still present in samples CI and C2 but is shifted of approximately 
28 cm^^ and 45 cm~^, respectively, upon compression. The ring-like mode associated to the threefold ring remains 
therefore Raman-active upon larger compression than the ones associated to the fourfold rings. This is certainly due 
to the fact that the Si-O-Si angles in the threefold rings have less possibility to change than in the fourfold rings. 
Indeed we checked that the Si-O-Si angles of the threefold ring in sample C, CI and C2 change by less than 3° 
upon compression, whereas changes of more than 20° were noticed in the case of the Si-O-Si angles of the fourfold rings. 



V. CONCLUSION 

We have presented results for VV and VH Raman spectra for uncompressed and compressed amorphous silica 
models generated by ab initio molecular-dynamics simulations. To model scattering in these materials we have 
considered the bond-polarizability mechanism. This mechanism gives rise to VV and VH spectra in good agreement 
with experimental data in shapes and peak positions. Analysis of the VV partial Raman spectra confirms that the 
Di and D2 peaks are due to the presence of four and threefold rings respectively, as suggested by Galeener. The 
association of ring-like modes to the peaks was done by using two approaches, one based on the atomic motions of 
the atoms belonging to the ring and the other one based on the calculation of the partial Raman spectra for atoms 
belonging to the ring. Both approaches give qualitatively the same results and show that not all rings necessarily 
produce a strong Raman activity. 

The Raman spectra obtained for compressed samples were also computed and are in very good agreement with 
experimental results by Hemley et al.^'^ for densified samples. As long as the tetrahedral network is maintained, we 
find that the Raman peak around 500 cm^^ is not necessarily due to modes associated to the fourfold rings although 
these rings are still present in the glass sample. However the D2 line still exists upon compression and can still be 
associated to the threefold rings. 
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These results suggest that (i) there exist coUectivc ring-Uke modes (involving several rings), (ii) the presence of the 
Di and D2 lines can be associated to the presence of fourfold and threefold rings in the glass sample, (iii) the absence 
of the Di and D2 lines does not mean that there are no fourfold and threefold rings in the sample, (iv) the area under 
the Di and the D2 peaks are not simply proportional to the concentration of small rings in the Si02 sample, as it has 
been previously proposed. 

As a conclusion, we were able to show that a combination of classical models (the molecular-dynamics simula- 
tions for the generation of the amorphous samples and the scattering models) and of first-principles calculations 
(the ab initio computation of the dynamical matrices) can be used to successfully compute the Raman spectra of 
amorphous silica and to extract detailed information about the Raman signature of structural units such as small rings. 
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TABLE I: Frequency, Raman intensity and localisation entropy for the most Raman-active modes in the 470-515 cm ^ frequency 
range (Di line) and in the 550-730 cm~^ frequency range (D2 line) for the three Si02 samples A, B and C. 

Di line 



sample A 


sample B 


sample C 


frequency 


Raman 


localisation 


frequency 


Raman 


localisation 


frequency 


Raman 


localisation 


(cm^"*") 


intensity 


entropy 


(cm-i) 


intensity 


entropy 


(cm^^) 


intensity 


entropy 




(a.u.) 






(a.u.) 






(a.u.) 




507 


5.76-10"'' 


0.801 


498 


5.32-10"-' 


0.806 


406 


7.16-10"'' 


0.831 


500 


4.86-10"^ 


0.804 


474 


4.38- 10"'^ 


0.678 


488 


6.76-10"=' 


0.777 


474 


4.42-10"^ 


0.757 


512 


4.37-10"^ 


0.811 


486 


5.36-10-^ 


0.804 








489 


2.09-10"^ 


0.836 


520 


3.96-10"^ 


0.823 








521 


1.87-10"=' 


0.847 


501 


2.20-10"^ 


0.872 








486 


1.77-10"^ 


0.847 








D2 lino 


sample A 


sample B 


sample C 


frequency 


Raman 


localisation 


frequency 


Raman 


localisation 


frequency 


Raman 


localisation 


(cm^"') 


intensity 


entropy 


(em^''") 


intensity 


entropy 


(em-i) 


intensity 


entropy 




(a.u.) 






(a.u.)' 






(a.u.)' 




634 


7.61-10"' 


0.81 


602 


1.05-10"^ 


0.82 


619 


1.53-10"^ 


0.75 


661 


2.88-10"-^ 


0.82 


609 


7.28-10"=' 


0.73 


652 


2.82-10"^ 


0.85 


606 


1.68-10"-^ 


0.82 


639 


6.9M0"=' 


0.72 


665 


2.12-10"^ 


0.79 


639 


1.64-10"^ 


0.85 


596 


5.51-10-^ 


0.83 









TABLE II: Pfting and Raman ratio (see text for definitions) for the most Raman-active modes in the 470-515 cm ^ frequency 
range for the three Si02 samples A, B and C. 

Sample A 



Frequency 




Pn 


ing 






Raman ratio 




(cm-i) 


ring 1 


ring 2 


ring 3 


ring 4 


ring 1 


ring 2 


ring 3 


ring 4 


507 


3.16 


2.35 


0.76 


1.35 


0.83 


0.70 


-0.32 


0.42 


500 


1.11 


2.60 


2.63 


0.40 


0.14 


0.59 


0.76 


0.21 


474 


0.30 


0.41 


0.44 


1.96 


0.08 


0.15 


0.16 


0.33 



Sample B 



Frequency 




Pning 




Raman ratio 


(cm-i) 


ring 1 


ring 2 


ring 3 


ring 1 


ring 2 


ring 3 


498 


0.39 


0.63 


3.81 


0.04 


0.23 


0.89 


474 


5.24 


0.35 


0.16 


0.95 


-0.04 


-5 10"^ 


513 


0.47 


2.99 


0.85 


0.19 


0.89 


-0.43 


489 


0.61 


1.23 


1.31 


0.11 


-0.52 


0.79 



Sample C 



Frequency 






Pr 


ing 










Raman ratio 






(cm""'") 


ring 1 


ring 2 


ring 3 


ring 4 


ring 5 


ring 6 


ring 1 


ring 2 


ring 3 


ring 4 


ring 5 


ring 6 


496 


0.50 


0.80 


0.98 


1.12 


2.62 


0.53 


0.22 


0.09 


-0.08 


0.18 


0.49 


0.23 


488 


3.99 


0.18 


0.76 


0.67 


0.53 


3.93 


0.91 


0.18 


-0.33 


-0.29 


0.05 


0.89 


486 


0.79 


1.59 


2.32 


2.55 


1.55 


0.61 


0.46 


0.66 


0.68 


0.73 


-0.52 


0.43 


520 


0.67 


1.37 


2.88 


1.17 


2.48 


0.32 


0.30 


-0.83 


0.99 


0.56 


0.90 


0.23 
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TABLE III: Number of common atoms in the fourfold rings of samples A, B, C, CI and C2. 





ring 


Sample A 
1 ring 2 ring 3 


ring 4 


ring 


Sample 
1 ring 2 


B 

ring 3 


ring 


Samples 
1 ring 2 ring 


C, CI and C2 
3 ring 4 ring 5 


ring 6 


ring 


1 




3 










1 



















3 


ring 


2 


3 










1 
























ring 


3 





























3 


3 





ring 


4 

























3 




1 





ring 


5 






















3 


1 







ring 


6 
















3 

















TABLE IV: PRing and Raman ratio (see text for definitions) for the most Raman-active modes in the 550-730 cm ^ frequency 
range for the three Si02 samples A, B and C. 



Sample A 



Frequency 


-PRing 


Raman ratio 


(cm-i) 


ring 1 


ring 1 


634 


3.33 


0.83 


661 


3.16 


0.97 


606 


0.69 


0.24 


639 


1.71 


0.83 



Sample B 



FrcqTiciic:v 


FRins 


Raman ratio 


(cm 


ring 1 


ring 2 


ring 1 


ring 2 


602 


1.47 


1.74 


0.29 


0.55 


609 


0.98 


6.24 


-0.09 


0.98 


639 


5.52 


1.43 


0.99 


-0.62 


596 


0.62 


1.08 


0.42 


0.38 



Sample C 



Frequency 
(cm-i) 


-PRing 

ring 1 


Raman ratio 
ring 1 


619 


5.65 


0.88 


652 


1.49 


0.78 


665 


0.58 


0.54 



TABLE V: Upper table: Frequency, Raman intensity and localisation entropy of the most Raman-active modes in the 490-540 
cm~^ frequency range for sample CI and in the 505-585 cm~^ frequency range for sample C2. Lower table: Frequency, Raman 
intensity and localisation entropy of the most Raman-active modes in the 605-770 cm~^ frequency range for samples CI and 
C2. 

Di lino 



Sample CI 


Sample C2 


Frequency 


Raman 


Localisation 


Frequency 


Raman 


Localisation 


(cm-i) 


intensity (a.u.) 


entropy 


(cm-i) 


intensity (a.u.) 


entropy 


512 


3.57-10-^ 


0.86 


546 


3.42-10"^ 


0.85 


503 


6.52-10"^ 


0.81 


534 


2.15-10"^ 


0.86 






Da 


line 






Sample CI 


Sample C2 


Frequency 


Raman 


Localisation 


Frequency 


Raman 


Localisation 


(cm-i) 


intensity (a.u.) 


entropy 


(cm-i) 


intensity (a.u.) 


entropy 


646 


1.54-10"^ 


0.79 


664 


1.48-10"^ 


0.79 


677 


4.56-10"^ 


0.83 


698 


9.65-10"'^ 


0.83 


632 


2.42-10"^ 


0.83 


643 


5.89-10"'^ 


0.87 








608 


5.80-10"^ 


0.80 
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TABLE VI: Pmng and Raman ratio (see text for definitions) for the most Raman-active modes in the 490-540 cm ^ frequency 
range for sample CI (upper table) and in the 505-585 cm"^ frequency range for sample C2 (lower table). 

Sample CI 



Frequency 
















Raman ratio 






(cm-i) 


ring 1 


ring 2 


ring 3 


ring 4 


ring 5 


ring 6 


ring 1 


ring 2 


ring 3 


ring 4 


ring 5 


ring 6 


512 


1.15 


1.15 


1.43 


1.28 


2.22 


0.58 


0.30 


0.26 


0.30 


0.18 


0.27 


0.21 


503 


3.35 


0.32 


1.02 


0.86 


0.78 


3.17 


0.82 


0.22 


-0.39 


-0.32 


0.16 


0.81 



Sample C2 



Frequency 

(cm^"'") 


ring 1 


ring 2 


Pr 
ring 3 


ing 

ring 4 


ring 5 


ring 6 


ring 1 


ring 2 


Raman 
ring 3 


ratio 
ring 4 


ring 5 


ring 6 


546 
534 


0.78 
1.85 


0.50 
0.24 


0.88 
0.85 


1.22 
1.39 


0.73 
1.09 


1.71 
1.00 


0.20 
0.40 


0.17 
-3 10-3 


0.19 
0.22 


0.24 
0.13 


0.03 
0.23 


0.34 
0.24 



TABLE VIL ftiing and Raman ratio (see text for definitions) for the most Raman-active modes in the 605-770 frequency range 
for samples CI (upper table) and C2 (lower table). 



Sample CI 



Frequency 


-PRing 


Raman ratio 


(cm-i) 


ring 1 


ring 1 


646 


4.83 


0.85 


677 


1.13 


0.48 


632 


1.18 


0.72 


Sample 


C2 


Frequency 


-Pfting 


Raman ratio 


(cm-^) 


ring 1 


ring 1 


664 


3.53 


0.79 


698 


0.65 


-0.21 


643 


0.98 


0.41 


608 


1.13 


0.51 
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FIG. 1: Frequency dependence of the effective density of states G{uj) obtained from ab initio (solid line) simulations and 
compared to neutron scattering experiments from Refi— (filled circles). The calculated G{lo) has been obtained using the 
experimental values for the neutron scattering length factors 6si=4.149 fm and &o=5.803 fm and a Gaussian broadening of 
width 2a— 35 cm~^. 
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FIG. 2: (a): Frequency dependence of the VV Raman intensity calculated using the ah initio eigen modes and the bond 
polarizability model (solid line). The dashed line is the experimental spectrum from Refp-. (b): Frequency dependence of the 
VH Raman intensity calculated using the ah initio eigen modes and the bond polarizability model (solid line). The dashed line 
is the experimental spectrum from Refi-. The calculated and experimental spectra have been normalized to the same total 
area. In the calculated spectra, the line shape is assumed to be Lorentzian and the line width is fixed to 25 cm~^. 




FIG. 3: Distribution of the rings size in the three amorphous silica samples A, B and C. 
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Sample C - mode 488 cm 
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FIG. 4: Ivv- Total VV Raman spectra of sample C (bold line). Imodc: VV Raman spectra in which the 488 cm~^ mode of 
sample C has been removed (thin line). laing: VV Raman spectra in which the eigenmode components of atoms belonging to 
ring 1 (upper graph), to ring 3 (center graph) or to ring 6 (lower graph) for the 488 cm~^ mode of sample C have been set to 
zero (dashed lines). 
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Sample B - mode 609 cm 
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FIG. 5: IvV- Total VV Raman spectra of sample B (bold line). I,nodc: VV Raman spectra in which the 609 cm^^ mode of 
sample B has been removed (thin line). laing: VV Raman spectra in which the eigenmode component of atoms belonging to 
ring 2 (upper graph), or to ring 1 (lower graph) for the 609 cm~^ mode of sample B have been set to zero (dashed lines). 
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Angles [degrees] 



FIG. 6: Distribution of the Si-O-Si angles (a) and of the 0-Si-O angles (b) of sample C at different densities: C at 2.20 g/cm^ 
(bold line), CI at 2.45 g/cm^ (dashed hne), C2 at 2.67 g/cm^ (dotted hne) and C3 at 4.28 g/cm^ (solid Une). 
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FIG. 7: Vibrational density of states of the sample C at different densities: C at 2.20 g/cw? (bold line), CI at 2.45 g/cm^ 
(dashed line), C2 at 2.67 g.cm^ (dashed-dotted line) and C3 at 4.28 g/cvc? (solid line). A Gaussian broadening of 2cr = 35 
cm~^ has been used. 
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FIG. 8: Frequency dependence of the normalized VV Raman reduced intensity (left panel) and of the VH Raman intensity 
(right panel) for the sample C at different densities: C at 2.20 g/cm^, 01 at 2.45 g/cm^, 02 at 2.67 g/cm^ and 03 at 4.28 
g/cm^. The spectra have been shifted vertically for clarity. 



